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Abstract 

The detectability of gravitational wave (GW) radiation from accretion disks is discussed based on 
various astrophysical contexts. In order to emit GW radiation, the disk shape should lose axial symmetry. 
We point out that a significant deformation is plausible in non-radiative hot accretion disks because of 
enhanced magnetic activity, whereas it is unlikely for standard- type cool disks. We have analyzed the 3D 
magnetohydrodynamical (MHD) simulation data of magnetized accretion flow, finding non-axisymmetric 
density patterns. The corresponding ellipticity is e ~ 0.01. The expected time variations of GW radiation 
are overall chaotic, but there is a hint of quasi-periodicity. GW radiation has no interesting consequence, 
however, in the case of close binaries, because of very tiny disk masses. GW radiation is not significant, 
cither, for AGN because of very slow rotation velocities. The most promising case can be found in gamma- 
ray bursts or supernovae, in which a massive torus (or disk) with a solar mass or so may be formed around 
a stellar-mass compact object as the result of a merger of compact objects, or by the fallback of exploded 
material towards the center in a supernova. Although much more intense GW radiation is expected before 
the formation of the torus, the detection of GW radiation in the subsequent accretion phase is of great 

importance, since it will provide a good probe to investigating their central engines. 
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1. Introduction 
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Gravitational wave (GW) radiation is expected to be detected within 10 years or so, thereby opening a new window 
to probe extremely high-density objects and the early universe (see, e.g. Thorne 1987, 1997; Blair 1991; Tsubono 
et al. 1997). When detected, the impact will be enormous and the understanding of astrophysical objects should 
inevitably undergo revolutionary development. It is thus of great importance to discuss at this moment what kinds 
of objects can emit GW radiation. In the present study, we consider if it is feasible to detect GW radiation from 
accretion disks in general contexts. 

At first glance, it is unreasonable to expect GW radiation from accretion disks, since, while the generation of 
GW requires the presence of quadruple moment, the usual disk models are constructed under the assumption of 
axisymmetry. Actually, the axisymmetric assumption is thought to be quite generally satisfied (see, e.g., Kato et al. 
1998 for a review of various disk models). Then, how and under what circumstances can non-axisymmetry arise? 

A non-axisymmetric disk structure is often discussed in the context of accretion disks in close binary systems. 
Since the shape of the Roche lobe is not totally spherically symmetric, we expect a disk there to lose axial symmetry, 
especially when the disk size is large, close to the size of the Roche lobe (Paczyhski 1977). Another good indication of 
non-axisymmetric structure is frequently discussed based on hydrodynamical simulations; e.g. Sawada et al. (1986). 
Note that the presence of spiral patterns on a disk was discovered through the observing technique of Dopplcr 
tomography (Steeghs et al. 1997). As we will see later, however, because disks in close binary systems have tiny 
mass, the production of strong GW radiation is unlikely. 

Alternatively, strong magnetic fields, if they ever exist, could greatly modify the disk structure. This is the 
subject of the present study. As was first discussed by Shakura and Sunyaev (1973) and later demonstrated by many 
authors through MHD simulations (Matsumoto 1999; Stone et al. 1999; Machida et al. 2000; Hawley et al. 2001), 
the magnetic energy can be amplified by the number of MHD instabilities together with differential rotation up to 
the value of p ma , g /p ga , s — 0.01-1 (e.g. Machida et al. 2000). The corresponding viscosity parameter is a ~ 0.01 — 0.1. 
This is just in the range that we require to account for the observations of dwarf-nova outbursts (Cannizzo 1993) 
and X-ray nova eruptions (Mineshige, Wheeler 1989). 

In standard-type disks, however, it is hard to believe a large influence of magnetic fields on the shape of the disk, 
since the magnetic energy is, at most, comparable to the internal energy of the gas, which is much less than the 
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kinetic (rotational) energy of the gas as a consequence of efficient radiative cooling of the disk material. In other 

words, the magnetic pressure cannot overcome the gravitational force in such radiation-dominated accretion flow. 

In advection-dominated regimes (see Kato et al. 1998 for a review), in contrast, the internal energy of gas can be 

comparable to its kinetic energy and potential energy because of a low radiation efficiency. We can then expect a 

large influence of magnetic fields on the dynamics of accretion disks. 

Another argument to support this idea comes from the observed complex variability commonly observed in many 
black-hole candidates during their hard state. In that state, the spectra are hard, and a hot accretion flow model 
(ADAF model) can fit the observations (Ichimaru 1977; Narayan et al. 1996). Although its origin is not yet 
established, many authors suggest that the variability could be caused by a sporadic release of magnetic energy 
triggered by magnetic reconnection and flares (Takahara 1979; Galeev et al. 1979). Without large magnetic field 
energy, it is difficult to account for substantial variations. This idea is consistent with the absence of large fluctuations 
during the soft state, when a standard-type disk seems to be present. It then follows that a non-axisymmetric disk 
structure may be generated by magnetic fields in cooling-inefficient regimes (i.e. adiabatic regimes) of disk accretion. 
In fact, the global MHD simulation of a disk exhibits a spatial inhomogeneous structure (Kawaguchi et al. 2000). 
Finally, note that large fluctuations are also observed in the high-luminosity state. Magnetic activities also seem to 
be enhanced in such a state (e.g. Mincshigc et al. 2000). 

We, here, consider GW radiation from accretion disks, the shapes of which are possibly influenced by large 
magnetic fields or other effects. In section 2 wc calculate the moments of inertia by using the MHD simulations data 
and see to what extent a deviation from an axisymmetric disk can be expected. We then discuss the detectability 
of GW radiation in various astrophysical contexts. The final section is devoted to discussion. 

2. Inhomogeneous Density Structure Created by Magnetic Fields 

Here, we analyze 3D MHD simulation data newly calculated by Machida and Matsumoto (2002). They started 
the simulation with a torus (Okada et al. 1989) threaded by weak toroidal fields (see Machida et al. 2000) in a 
pseudo-Newtonian potential (see similar calculations by Hawley, Krolik 2001). The size of the calculation box was 
taken to be 100 r g (with r g being the Schwarzschild radius) and the mass within the last stable circular orbit at 3 r g 
was removed as an inner boundary condition. We have confirmed that magnetic fields are amplified up to the value 
of Pmag/fgas — 0.1 within ~ 50 rotation periods at the reference radius of r = 50 r g (sec Machida, Matsumoto 2002 
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Fig. 1.. Snapshot of the density contours of a magnetized disk on the equatorial plane [at a calculation time of t = 32070(r g /c), see 
Machida, Matsumoto 2002]. This figure covers the region of 20 r g X 20 r g . The spacing between each contour is Alogp = 0.05. 

for more details). 

Figure 1 represents a snapshot of the density contours on the disk equatorial plane. A sort of spiral pattern (with 
the azimuthal number of the spiral being m = 1) is clear in this plot. As time goes on, this pattern shows a rotation 
around the center, thus producing a non-zero ellipticity. However, the pattern is not the same in the rotational 
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frame due to the presence of differential rotation and magnetic energy dissipation by magnetic reconnection. 

To evaluate the strength of GW radiation, we next calculate the moments of inertia using the MHD simulation 

data. Here, we use cylindrical coordinates, (r,tp,z), with the z-axis being taken as the rotational axis. Then, the 

moments of inertia are written as 

I xx = J p(r,<p,z)x 2 rdrdtpdz (1) 



and 



= J p(r, p, z)y 2 rdrd(pdz (2) 



Iyy — 

with x = r cos ip and y = ramp. Since the spiral density pattern rotates with time, such a magnetized disk can 
certainly produce GW radiation, which is rather chaotic, as is shown in figure 2. At the same time, we also notice 
a hint of quasi-periodic variations on timescales of <~ 200 r g /c in the variations of the cllipticity (e) and moments of 
inertia (I xx and I yy ). 

The different behavior of I xx and I yy yields an ellipticity, which is calculated by 

_ Ixx — Iyy /„n 
6 ~ T + J ' ^ 

±xx r Iyy 

As is well known, non-axisymmetric disks (with e ^ 0) can emit GW radiation. In the next section we estimate 
how much energy can be extracted by the GW radiation and how large an amplitude of GW radiation is expected 
for a simple case of a non-axisymmetric disk with an ellipticity of e <~ 0.01. It might be noted here that, precisely, 
although we need to sum up all of the contributions from each mass particle moving around the center, for the 
purpose of order-of-magnitude discussion presented in this paper, it will suffice to use the formula for an elliptical 
mass distribution. 

3. Possible GW Radiation from Accretion Disks 

3.1. Basic Formulae 

In the previous section, we have seen how a non-axisymmetric disk structure is created by global magnetic fields. 
We expect about a few percent of inhomogeneity in the magnetized disk. Wc then proceed the discussion and 
consider what the most ideal targets would be. 
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Fig. 2.. Time variations of the cllipticity (upper panel), and those of moments of inertia, I xx (by the solid lines) and I yy (by the 
dashed lines), respectively (lower panel). The units of the time and moment of inertia are r g /c (si 3 X 10~ 5 s for the central mass 
of M» = 3 Mq) and (M disk /M Q )(r d i sk /10 7 cm) 2 , respectively. We find e ~ 0.01. 

The energy loss rate due to the GW radiation is expressed by 



• 32 G 2 26 



(4) 
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whereas the expected amplitudes are on the order of 

h+ = - (h) (cos 2 9 + 1) cos(2ffi - 2<p ) (5) 

and 

h x = -2(h) cos 9sm(2Qt- 2ip ) (6) 
with a coefficient of 

(h) = {2eGm 2 )/(c 4 d), (7) 

and ipo being a constant. Note that the effective amplitude is h c s — VN (h) with N being the number of the cycles 
with regard to the quasi-periodicity during the observing time. 

3.2. Case of Binary Accretion Disks 

For disks around neutron stars or black holes in close binary systems, typical parameters are a disk mass of 
-Mdisk ~ 1O~ 8 M0, a disk size of rdisk ~ 10 10 cm, and an orbital period at the disk rim of 



where m* = M*/M Q with M* being the mass of the central object, and no = rdi s k/10 10 cm. We then find 
where mdisk = A^disk/-^© and e_2 = e/10 -2 . The timescale of the energy loss from the disk would be 

-E-rot 



l-Egwl 



Since the characteristic timescale greatly exceeds the age of the universe, GW radiation is totally unimportant for 
the evolution of binary disks. The expected GW amplitude is 

for TV = 1000. This is far below the current detection limit. There is no hope to detect GW radiation from the 
binary accretion disks. 
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3.3. Case of Massive Disks in GRBs or HNe 

Next, we consider a massive disk around a stellar-mass compact object. It has been proposed that a massive torus 
or disk (with solar mass or so) is possibly formed around a stellar-mass compact object after the merger of compact 
objects, a promising possibility to produce gamma-ray bursts (e.g., Narayan et al. 1992; Kohri, Mineshige 2002). 
Likewise, fallback of exploded material towards the center may also form a massive torus in a supernova (SN) or 
hypernovae (HN) (Wooslcy 1993). For recent reviews of GRBs, see Piran (2000); Meszaros (2001), and references 
therein. We thus expect much more intense GW radiation in these cases than in the case of close binaries. 

The typical parameters are ~ M and r^sk ~ 10 7 cm; and the orbital period at the disk rim is 

(f) _1 V[s], (12) 

with r'j = rdi s k/10 7 cm. That is, the modulation frequency is / = 2/P ~ 200 Hz or so. From the formula we find 
the energy loss rate to be 

4w = -4.5 x io«!^a (I^) 3 [erg s" 1 ]. (13) 



The timescale of the energy loss from the disk is estimated to 

i gw ^10 4 -4^(^)" 2 [s]. (14) 

Since the lifetime of a disk in one GRB is several tens to thousands of seconds, at longest, GW radiation is not 
very important. Note, however, that more deformation of the disk shape is plausible if the massive torus is unstable 
against self-gravity. Then, GW radiation may affect the secular evolution of the disk. The effective amplitude would 
be 

-22 /irr m diske-2 (m t 



M ,1.3x10-^3=^^; , (15) 

where N 3 = TV/1000 and we took N <~ 10 s / 0.01 s = 10 3 . This is marginally detectable by the LIGO, as long as a 
GRB occurs within a distance of <~ 10 Mpc. Hence, there is a possibility of detecting GW radiation from hypernovae 
or GRBs, although much more enhanced GW radiation is expected in a violent phase (i.e. SN explosion or tidal 
destruction of a star by a black hole before a merging) preceding the formation of a massive disk (e.g. Yamada, Sato 
1995; Shibata, Uryu 2000). Nevertheless, the detection of GW radiation from the disk accretion phase is crucial, 
since it would continue for the burst durations, thus providing an essential clue to understanding what is actually 
going on at the center of SNe and GRBs. 
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3.4- Case of AGN Accretion Disks 

Finally, let us consider accretion disks in AGN. There is a convenient scaling law describing the mass distribution 

and time scales of AGN accretion disks. Since the surface density does not strongly depend on the black-hole mass, 

as long as the radial distance to the center is expressed in terms of the Schwarzschild radius (r g oc Mbh) and the 

mass- accretion rate is scaled by the Eddington rate (^Edd/c 2 with LEdd being the Eddington luminosity), we have 

the following approximate scaling laws: 

A/disk ~ 7rrL k £ « A*bh (16) 

and 

I ~ M disk r 2 disk « M* H . (17) 

These scalings exactly hold, if the disk is radiation pressure-dominated and if the disk size is proportional to the 
black-hole mass. Numerically, we thus estimate 

M disk ~lO- 8 (^)V «lO 7 O'Me (18) 

and 

r disk - 10 10 (^) [cm] « 0.1 (^t)M. (19) 
Likewise, the orbital periods and angular frequencies can be scaled as 

P~10 10 (^p) [ S ]ocM BH (20) 
(with tobh = Mbk/Mq) and, hence, Q cx Therefore, we find 

( sw a/ll 4 aMj| H , (21) 

that is, the characteristic timescale is mass-independent and thus again exceeds the age of the universe [see equation 
(|lo|)1. (The timescale may depend on the black-hole mass, if the disk size is not proportional to the black-hole mass; 
but, even for such a case it is highly unlikely that GW radiation significantly affects the time evolution of AGN 
disks.) 

How about the effective amplitudes? By using equation (nfi), we easily obtain 



h+ cx m 2 /d(x M| H , (22) 
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leading to the estimation 

U = U .^(=)'( 5 i ; )-'. (23, 

Thus, the amplitude is not small, as expected, because of a massive disk. However, the frequency is only / = 2/P ~ 
2 x 10~ 10 Hz; that is, there are no GW detectors planned which are capable of detecting such extremely low GW 
radiation. 

Because the disk is rotating more rapidly in the inner portions, there may be a chance of detecting GW radiation 
from the innermost part of the AGN disks. Let us consider the innermost part with a radius of r . Then, the mass and 
the moment of inertial are reduced by factors of (r ou tAo) 2 and (rdi s k/?"o) 4 , respectively, while the angular frequency is 
increased by a factor of (^disk/m) 3 / 2 , as long as the Kepler rotation law is adopted. By setting r ~ 10 r g <~ 10~ 3 rdi s k, 
we obtain 

and the frequency is 

/~6xl0- 6 (^) -1 [Hz]. (25) 
Still, it is difficult to detect GW radiation for AGN located at d ~ 10 Mpc by LISA. 

4. Discussion and Conclusion 

We have discussed the detectability of GW radiation from accretion disks based on the assumption that disks are 
magnetically dominated, and thus their structure is inhomogeneous, leading to the emergence of non-axial symmetry. 
In the case of close binaries, however, GW radiation is totally negligible because of their very tiny disk masses. In 
the case of AGN, GW radiation is not significant, cither, because of slow rotation velocities. Instead, we may 
concentrate on the innermost part of the AGN disks, which has much faster rotation velocities. In fact, the inner 
hot part of AGN disks seems to be in a non-radiative state, thus producing hard X-ray emission. Yet, the detection 
of GW radiation is not very feasible with the currently planned GW detectors, unless additional mechanisms work 
to enhance the disk deformation. 

The most probable case may be found in gamma-ray bursts and supernovae (or hypernovae), in which a massive 
disk (with solar mass or so) is suggested to be formed around a stellar-mass compact component. Note, however, 
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that the time variations of GW radiation are not predictable; they are highly chaotic, although there is a hint of 
quasi-periodicitics. Also note that much more intense GW radiation is expected before the formation of a massive 
torus. Nevertheless, the detection of GW radiation from such a massive disk is of great importance to probe the 
central engines of the GRBs and SNe (HNe). 

The accretion disks can emit electro-magnetic (EM) wave radiation as well as GW radiation. The former might 
be even more important than the latter, when considering the long-term evolution of the disk, since whereas the 
emissivity of GW radiation is proportional to f2 4 , that of the EM wave radiation is proportional to SI 2 . One might 
think that the disk rotation may then be efficiently decelerated by the EM wave radiation so that significant GW 
radiation cannot be expected. However, this does never occur, since the removal of angular momentum of a disk 
gas results in faster rotation, rather than slower rotation, unlike the case of rotating stars. This is because when 
the disk material loses angular momentum, £ (oc \/r), the disk gas moves inward because of the reduced centrifugal 
force (which is oc I 2 ), leading to more rapid rotation with higher angular frequency (note that S} oc r~ 3 / 2 ). This 
makes a good contrast with the case of rotating neutron stars, in which the removal of angular momentum results 
in a spin-down of the star. 

Further, we find it unlikely that the evolution of the disk is totally modified by the EM wave radiation, since the 
magnetic energy, which produces EM wave radiation, is only ~ 1% of the gravitational energy (which is equal to 
the rotational energy by the Virial theorem). To be more precise, the Alfven wave is more important than EM wave 
radiation in terms of the energy loss from the disk. The Poynting flux due to Alfven waves propagating along the 
poloidal magnetic field is 



where v is on the order of the rotation velocity of the footpoint of the poloidal field anchored by the disk; that is 
v = rCl. On the other hand, the energy loss due to the EM wave radiation is at most 



Eah ~ B 2 vr 2 , 



(26) 



E EM ~ (Sl 2 Br 3 ) 2 /c 3 , 



(27) 



since the magnetic moment is of the order of Br 3 at maximum. Thus, we have 



Eem ~ E AU (rSl/c) 3 < E AU . 



(28) 
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The characteristic time of the disk evolution due to the emission of Alfven wave is then 



*Aif ~ {GM.p/rjB- 2 ^ 1 > ft 



-l 



(29) 



[recall B 2 /(GM*p/r) ~ 0.01]. Note that this seems to be highly underestimated, since the magnetic fields have 
bisymmctric spiral patterns (Machida, Matsumoto 2002) so that the magnetic moment should be much less than 
what we assumed above. We can thus safely conclude that the EM effects are not very significant, and, thus, our 
argument concerning the GW radiation from the disk is not basically altered. 

Since different parts of the disk material rotate with different angular frequencies (ft), the actual GW radiation 
could be the superpositions of many modes which vary on a variety of periods. This leads to a reduction in the 
effective GW amplitudes estimated in the previous section. Whether this is the case or not strongly depends on the 
magnetic-field distribution in the disk. (Such a problem does not occur in the GRB case, since only a compact disk 
has been considered there.) If magnetic fields are dominant in the innermost region (within, say, several hundreds of 
Schwarzschild radii), clear spiral patterns are only appreciable in a relatively narrow region. Unfortunately, however, 
it is still difficult to calculate the long-term evolution of the magnetized disk with sufficiently large spatial dimension. 
Further, it is currently impossible to simulate two-phased disk corona structure with proper treatments of magnetic 
fields and radiation transfer. We need to await until it becomes possible to calculate the coupling between the 
radiation, magnetic fields, and matter in accretion flow. 

Strong GW radiation from the gamma-ray burst was discussed by van Putten (2001) in a slightly different context. 
He conjectured that magnetic fields with poloidal topology equivalent to pulsar magnetospheres provide tight coupling 
between a Kerr hole and suspended accretion torus, thereby radiating GW by extracting the spin-energy of the black 
hole in the presence of non-axisymmetry. In the present study, in contrast, we do not assume the presence of rapidly 
spinning black holes, nor tight coupling between the black hole and the torus, which makes a distinction between 
his model and ours. Our model should work in any disks with hot gas, in general, and does not require any special 
condition. 

We wish to finally note that for a massive disk whose mass is comparable to, or exceeds, that of the central object, 
the torus itself becomes gravitationally unstable, thus giving rise to non-axisymmetry, even without magnetic fields 
(e.g. Woodward et al. 1994). Then, the corresponding e value will be much increased: say, e ~ 1 in the extreme 
case of fragmentation of the disk matter into several pieces. This may be relevant in the case of a compact massive 
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disk as in GRBs and HNe and strengthens our conclusion. 
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